ROTOR INTERACTION IN THE ANNULUS BILLIARD 



PETER BALINT AND SERGE TROUBETZKOY 

Abstract. Introducing the rotor interaction in the integrable sys- 
tem of the annulus billiard produces a variety of dynamical phe- 
nomena, from integrability to ergodicity. 



1. Introduction 

The rotor model (see (1.1) below) has been introduced in [MLL] and 
further studied from the physical point of view in [LLM]. The main 
motivation for replacing the usual elastic reflection in billiard models 
with such an interaction is that this way the a priori independent mo- 
tion of point particles can be coupled. Despite of the physical interest 
in the subject, there is only a very little amount of mathematically 
rigorous discussion. This is related to the difficulty of the problem - 
even the simplest possible hyperbolic candidates, planar Sinai billiards 
with circular scatterers and the related Lorentz gas models turned out 
to be rather complicated. In [BK] a detailed analysis has been given 
for a zero measure set of trajectories that possess special symmetries. 

In the present article we investigate the rotor interaction in the con- 
text of elliptic billiards. Consider the motion of point particles in the 
annulus billiard, i.e. in the planar region bounded by two concentric cir- 
cles. The one particle system is well known to be integrable ([T] p. 26), 
while the dynamics of several point particles is the direct product of 
the individual regular systems. However, if we replace the inner (and 
possibly the outer) circular billiard scatterer with a rotating object, 
regularity may be destroyed. Three different models are investigated 
in this article. In Section 2 we consider the one particle system with the 
inner circular scatterer rotating and the outer one kept to be a billiard. 
For this system integrability persists. In Section 3 both scatterers are 
assumed to rotate and this setting destroys integrability. For a positive 
measure set of parameters velocity evolution is minimal (Theorem 3.7) 
and the two dimensional system of our point particle is ergodic (The- 
orem 3.15). Finally, in Section 4 the motion of two point particles is 
discussed in case the inner scatterer is rotating and the outer one is kept 
to be a billiard. Though this problem is significantly more complicated, 
it is possible to find an open set of parameters and, correspondingly, a 
positive measure set of initial conditions, that this system behaves as 
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the one discussed in Section 3 for an arbitrary long finite time (Theo- 
rem 4.4 and Remark 4.6). Prom a physical point of view this result is 
probably the most interesting as it indicates, as claimed in the physics 
literature (see [MLL]), the following phenomena: the rotor interaction 
switches on non-trivial coupling among the particles that would follow 
independent (in our case integrable) patterns if they were studied in 
the original biUiard setting. 

Throughout the paper wc consider the motion of point particles in 
the annular region bounded by two concentric circles, the outer hav- 
ing radius one, the inner radius R < 1. Interaction with the circular 
scatterers is either an elastic reflection (as in billiard models) or the 
scatterer is assumed to rotate with some angular velocity to and interact 
with the point particle according to the laws 

(1.1) v't = vt - Y:^{vt - ru); 

2 , 

ru> — ru> -\ [vt — ru>). 

1 + r] 

Here the coordinates y„ and Vt are the normal and tangential com- 
ponents (with respect to the circle of coUision) of the particle velocity 

and the primes indicate post-collisional values. Furthermore r is 1 for 
the outer and R for the inner circular scatterer, and ri = where 
is the moment of inertia for the scatterer in question and m is the 
mass of the point particle. 

Remark 1.1. Note that for the same velocity vector v the splitting 
into tangential and normal components is different for the inner and the 
outer scatterer. In Sections 2 and 4 it is only the inner scatterer that is 
assumed to rotate, and thus the components Vn and Vt in these sections 
always refer to the splitting at the inner scatterer. In Section 3 both 
scatterers rotate and it is more convenient to use different notation, 
see the beginning of Section 3.1. 

Convention 1.2. The evolution of the coordinates in Formula (1.1) 
will be referred to as the base. 

In each section the dynamics of the base will be studied before turn- 
ing to the whole system, which is a skew product over the base, some- 
times referred to as fibers. 

2. One particle, one rotating scatterer 

This section is devoted to a model which is closest to the original 
integrable annulus billiard, namely the motion of one point particle 
with only the inner scatterer rotating. For certain initial conditions 
the point particle never reaches the inner circular scatterer. In this 
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invariant set motion is that of a point particle in a circle, a system 
well-known to be integrable. We consider the complement. 
At the inner scatterer we have 

(2.1) v\ = vt-^{vt-Ruj)- 

2 

Rjj' ^ Ru;-\ (vt - Ew) 

l + rj 

which is just the inner version of Equation (1.1), see Remark 1.1. 

Lemma 2.1. (1) If a trajectory collides once with the inner scatterer, 
it collides with it infinitely often. More precisely, there is an alternat- 
ing sequence of collisions, every first with the inner and every second 
with the outer scatterer. 

(2) There are three integrals of motion. 

(3) The velocity motion is periodic of period 2. 



Proof. A trajectory which starts from the inner circle is reflected elas- 
tically in the outer circle, thus it returns to the inner circle. Iterating 
this argument proves statement (1). 

Consider two consecutive reflections at the inner scatterer with one 
billiard reflection at the outer one in between. Let us denote by the 
triples v„(l), ft(l), a;(l) and Vn{'2),Vt{2),uj{2) the incoming velocities at 
the two consecutive inner reflections, and let us refer to the outgoing 
counterparts by primes. By the geometrical picture sketched above: 

^;„(2) = -w„(l)'; vt{2) = vt{l)' 

and furthermore a;(l)' = tf (2) as the inner scatterer does not interact 
while the particle collides outside. This, together with (2.1), gives, on 
the one hand 

Vr,{l) = i;„(2) = -Vr,{l)' = -Vr,{2)' 

which can be regarded as an integral of motion. On the other hand, 
to get Vt{2)' and uj{2)' we need to iterate the second and the third 
formulas in (2.1) twice. This gives 

u{2)' ^uj^l); vt{2)' ^vt{l) 

thus statement (3) is proved. 

Finally multiply the third equation in (2.1) with r] and add it to the 

second to get 

(2.2) Vt + T]Ruo = v[ + r]Ruj' = N. 

Furthermore, subtracting the third line from the second in (2.1) we get 

v't-Ruj' ^ -{vt - Ruj) 



4 



ROTOR INTERACTION IN THE ANNULUS BILLIARD 



which, together with (2.2) resuhs in 

(2.3) v1 + r]R^u'^ = v'^ + t^R^u'"^ = E. 

Thus we have calculated the other two integrals. □ 

Remark 2.2. Equation (2.2) is a multiple of the angular momentum, 
while (2.3) is the "tangential" part of the full kinetic energy (the normal 
part of the kinetic energy, i.e. the length of Vn is in itself conserved as 
the third integral). 

In the space of velocities {vt,uj), that can be regarded as R^, fixing the 
values (2.2) and (2.3) we get a line and a circle, respectively. These 
two curves intersect in ( at most) two points {vt, cu) and {v^, cu') . Velocity 
evolution is a period two cycle on these two points. 

Turning to the evolution in the fibers, we essentially get the behavior 
observed in the integrable annulus billiard. Let us parameterize the 
outer circle by the arclength parameter (f and denote the points of 
consecutive impacts on a trajectory by ipk,k G Z. Then 

^k+i = (Pk + P + P' 

where 

P — arctan— , P — arcsin 

Vn 

and the same for the primed quantities, see Figure 1. Thus the angles 
/3, P' G [— arccos(i?), arccos(-R)] depend smoothly and strictly monoton- 
ically on the ratios ^ and respectively. The consecutive (/^^-s follow 
an orbit of a circle rotation, which is irrational apart from a countable 
collection of one codimensional manifolds in the three dimensional pa- 
rameter space of integrals. To see this keep the values of E and N 
((2.3) and (2.2)) fixed and increase the value of f„ continuously, then 
P + P' decreases continuously and is rational only for a countable set 

of Vn-S. 

3. One particle, two rotating scatterers 

3.1. Coordinates, parameters and integrals of motion. In this 

section we will find that allowing both scatterers to rotate results in a 
behavior different from the integrable models discussed so far. 

As long as velocities are concerned we need to handle an a priori 
four dimensional system: the coordinates are, on the one hand, u>i,u!2 
for the angular velocities of the inner and outer scatterers, respectively, 
and, on the other hand, v. the two dimensional velocity vector of the 
point particle. Recall Remark 1.1: in contrast to the previous section, 
we will use notations Vn and Vt for the splitting of the velocity v into 
normal and tangential components at the outer scatterer. The sign of 
Vt has a physical meaning, it describes the direction of the "circular" 
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component of the motion of the point particle. On the contrary, the 
sign of Vn is in a certain sense irrelevant, the positivity /negativity of 
Vn merely tells us if we are just before/ after a collision with the outer 

circle, thus we fix Vn to be always positive. The symbols Vn, Vt arc used 
for the splitting of v into normal and tangential components at the inner 
scatterer. By simple geometrical observation they can be defined for 



vectors v which satisfy ^ < j^j^ as Vt = and Vn = y^n — ^-^-vj- 
Otherwise the trajectory starting from the outer scatterer reaches it 
again without interacting with the inner one. 

We use the notation ooi for the angular velocity and Oj for the mo- 
ment of inertia with i = 1,2 in the case of the inner and the outer 
scatterer, respectively. The rescaled moments of inertia are ?7i = 
and m = —■ In these coordinates, the interaction at the outer scatterer 
is given by: 

V, = Vt - ^^{vt - u;^); 
2 

(3.1) LO2 = UJ2 + — {Vt-UJ2)- 

By the convention on the sign of Vn made above, the normal velocity is 
preserved and straightforward calculation gives two further preserved 
quantities: 

(3.2) Vt + r]2UJ2 = v't + r]2UJ2, vj + 772^^2 = '^'t + 'n2^2 ■ 
As to interaction with the inner scatterer 

2?7i 

v't = Vt - -^^^^ {vt- ^1); 
2 

Rjjj\ = Ruji H (vt — Ru)i). 

1 + ?7i 
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Note, however, that Vt = Rvt and thus these equations are equivalent 
to 

Vt = '"t - ^ ^^ {vt - R ^i); 

2 

(3.3) R^ui = R^ui + (vt - R^uJi). 

1 + ??i 

Two preserved quantities for the inner collisions are 

(3.4) vt + r}iR'^uJi^v[ + rjiR^co[; + rjiR'^col ^ + rjiR^cof . 

As uji {UJ2) does not change at outer (inner) collisions, (3.2) and (3.4) 
imply the presence of two integrals of motion, the angular momentum 

(3.5) Vt + r]iR^u>i + 7)2(^2 — const. — N 
and the tangential energy 

(3.6) + rjiR'^Lof + r]2Ujl = const. = E{> 0). 

The term "tangential" refers to the fact that the full kinetic energy 

of the system is a third integral of motion, always greater than the 
tangential one. More precisely, dividing the full energy ^{mxP' + QiujI + 
62t<^2) by \'m we get 

(3.7) vl + v'l + r^iR^ul + 772^2 = const. = F(> E). 

Note that (3.2) ((3.4)) and (3.7) together imply that Vn {vn) is preserved 
during inner (outer) collisions. 

Remark 3.1. (1) One of our aims below is to understand how the 
dynamical behavior depends on the several parameters of the system. 
The constants rji and 772 will be referred to as physical parameters in 
contrast to the integrals of motion N, E and F. 

(2) To avoid complications we assume v ^ 0. Note that this as- 
sumption is reasonable in the following sense: v = is only possible 
for a codimension one submanifold of integrals of motion. To see this 
observe that if we had v = 0, (3.5) and (3.6) would fix the values of 
the other two coordinates ui and uj2 (up to two possibilities) and thus 
would create, via (3.7), a relation between the three integrals of motion. 

3.2. Evolution of velocities. For convenience we rescale our veloci- 
ties as 

(3.8) X := ^/r]^R^UJl, y := ^/rj^u;2, z := Vt, w := v^. 
In these rescaled quantities our integrals of motion are 

(3.9) x^ + y^ + z"^^ E; ^x + + z^N 
and 

(3.10) ^ + y^ + + = F. 

R^ 
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By (3.1) and (3.3) we can first consider how the velocity coordinates 
x,y,z evolve at collisions. With the convention on its sign the value of 
w{= Vn) is determined by (3.10) for each fixed value of the other three 
coordinates. In the rest of the section, unless otherwise stated, velocity 
evolution will always refer to the dynamics in the triple of x, y, z. 

Taking into account that x {y) does not change at outer (inner) 
collisions, in the rescaled coordinates the transformations (3.1) and 
(3.3) are reflections across the planes 

(3.11) y/fj^z - y = and ^/fjlz - x = 0, 

respectively. Furthermore, by (3.9) motion is restricted to the inter- 
section of a sphere and a plane, i.e. to a circle which will be referred 
to as the velocity circle. We will see that the base, see Convention 1.2, 
will consist of two copies of this circle. Note that these circles do not 
describe the physical position of the particle, on the other hand they 
give the state of the velocities. On the velocity circle the reflections 
across the planes (3.11) reduce to reflections across lines. For the outer 
collisions we have reflections across 

(3.12) /2(t)= f-^^t, v^t, + 
while for the inner ones across 



Let us denote the angle between these two lines by |. Depending on 
it we distinguish rational, irrational (and later on even Diophantine) 
cases. 

Lemma 3.2. Rationality is independent of the integrals of motion. On 

the other hand, there exists a countable set of codimension 1 subman- 
ifolds of physical parameters which lead to the rational case, otherwise 
(and thus for a full Lebesgue measure set of parameters) 7 is irrational. 

Proof. Straightforward calculation gives 

(3.13) cos| = ((l + r/r')(l + 772-^))-^. 

Thus the rationality of 7 depends only on the physical parameters. 
Consider the numbers cos(||) for positive integers p < q. Note that 
elements of this countable set A are all algebraic numbers. The base is 
rational iff the left hand side of (3.13) is equal to some number in A. 
This happens for a countable collection of one codimensional subman- 
ifolds of physical parameters. □ 

Convention 3.3. To agree with different standard conventions, in the 
geometric framework we will think o/S^ as [0, 27r), while in the number 
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theoretic framework we think o/S^ as [0, 1). In both cases we scale all 
E)^ valued quantities (s,(f,'y,a{s)) appropriately. This should cause no 
confusion. 

To give a description of the base dynamics we introduce the arclength 
parameter s as a coordinate on the velocity circle. For definiteness 

s = corresponds to one of the intersection points of the circle with 
the line (3.12). In this coordinate the outer and inner bounces, i.e. the 
transformations (3.1) and (3.3) are 

(3.14) s — > Si — —s (modi) and 

(3.15) s — > S2 = — s — 7 (mod 1), 

respectively. Take two copies, O and /, of the velocity circle describing 

the outgoing velocities just after the bounces at the inner and outer 
circles, respectively. The phase space of the velocity motion, or the 
base, is 

M = §^ X {1, 2} = / U O; M3s^{s,i) with s e §^ and i = 1, 2. 

The velocity evolution will be denoted by T : M ^ M. As an inner 
bounce is always followed by an outer one, for s G /, i.e. for (s, 1) the 
image Ts is always in O, more precisely T(s, 1) = (si, 2) (recall (3.14)). 
An outer bounce, however, may be followed by either an inner or an 
outer bounce, depending on the angle the point particle velocity makes 
with the normal vector of the outer scatterer. 

Lemma 3.4. (1) There is a (possibly empty) open set U C S"*^, such 
that for s & U (more precisely for s — (s, 2) with s & U), we have 
Ts = (si,2) (the image is the point defined by (3.14) in O), and for 
s E 0\U (more precisely for s = (s, 2) with s E E>^ \U) we have 
Ts = (s2, 1) (the image is the point defined by (3.15) in I). 

(2) The set U dS^ is invariant with respect to the reflection (3.15). 

(3) There is an open set of integrals of motion for which U is empty. 

Remcirk 3.5. We will refer to empty U as the alternating case, non- 
empty U as the non- alternating case. Both of them occur for open sets 
of integrals of motion, independently of the value of physical parame- 
ters. 

Proof of Lemma 3.4. 

(1) By Formulas (3.9) and (3.10) the arclength s on the velocity cir- 
cle determines z = Vt and w = Vn > continuously. The geometric 
condition for (s, 2) = s e O to reach the outer scatterer without touch- 
ing the inner one (and thus interacting with it) is ^ > j^^- Note 
that at tangential collisions with the inner scatterer is modified by 
(3.3), thus this corresponds to landing on /. 
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We define 

(3.16) e §M As) > Y^2^%') } • 

If s G O witli s eU (i.e. in case tlie next bounce is outside) tlien (3.14) 
applies, and thus Ts = (si, 2). Otlierwise (i.e. if s G \ (f/ x {2})) the 
next bounce is inside, thus the image is the point defined by (3.15) in 
/, or Ts — (s2, 1). The set ?7 C is the preimage of an open set by a 
continuous map, thus it is open. 

(2) To see that U is invariant under (3.15) note first that y docs not 
change during an inner collision. This, together with the presence of 
the integrals (3.9) and (3.10) implies that w'^ — ^-^-z^ is invariant for 
inner bounces, i.e. for the reflections (3.15). The statement follows as 
it is exactly the sign of this quantity that distinguishes the points of 
U. 

(3) The set U is empty if 

(3.17) z\s) < Y^yis) 

for all s. Consider the integrals of motion defined by (3.9) and (3.10) 
and denote A = ^ > 1. By (3.9) and z"^ are both less than or equal 
to E for all s. This, together with (3.10) gives 

> F-y'^-z^ -x^ - x^{R-'^ - 1) > F-E- (i?"^ - 1)E = 
= F - ER-^ = E(X - R-^) > (A - R-^)z^ 

for all s. This implies U is empty if A — > ^^P-, i-e. if A > {^ — 1), 
a condition satisfied for a set of parameters that has nonempty interior. 

□ 



Remark 3.6. (1) Actually U is empty iff we replace the strict inequal- 
ity in Formula (3.17) with a nonstrict one. The analysis of the equality 
case is substantially more complicated, and is a co-dimension one phe- 
nomenon in the set of parameters, thus we do not analyze it. 

(2) Note that the points s & U correspond to velocity configurations 
that cannot be realized at the inner circle. Thus for nonempty U the 
phase space is M d M: 

M^IUO^ ((§^ \U)x {1}) U X {2}) 

and the dynamics acts as 



T(0\(C/x{2})) = (§1 \ t/) X {1} = /; 
T(U X {2}) UTI = §^ X {2} = O; 
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with 





T{s, 


2) 


= (^2,1) 


for 


seS^\U; 


(3.18) 


Tis, 


2) 


= (^1,2) 


for 


s E U; and 




T{s, 


1) 


= (si,2) 


for 





which is well-defined and invertible by statement (2) of Lemma 3.4- As 
U has no inside copy sometimes we refer to U x {2} as U for brevity, 
this should cause no confusion. 

The main result on the base dynamics is the following theorem. 

Theorem 3.7. In the alternating irrational case the dynamics of the 
base is minimal, while in all the other cases velocity motion is always 
periodic. Moreover, all possible periods are even with length bounded 
above by a constant (which depends on the parameters). 

Proof. Consider the alternating case first. For s G O all odd iterates 
are in I and all even iterates are in O (and for s e / the other way 
round). Thus it is enough to consider the first return map onto O 
(which is exactly T^). By (3.14) and (3.15) this is a rotation of §^ by 
7, which is periodic/minimal exactly in the rational/irrational cases. 
(The odd iterates make the reflected image of the same rotation orbit 
in /.) 

Turning to the non- alternating case, i.e. to t/" 7^ 0, first note that 
U r\T~^U is an invariant set. To see this consider (s, 2) — s E Ur\T~^U. 
Then both Ts and T^s are in O. More precisely 

T{s, 2) = [si, 2) and T^{s, 2) = T(si, 2) = [s, 2), 

thus the set is invariant, furthermore, all points of it are periodic of 
period 2. 

Now take {s,2) ^ s e U \ T'^U. Then Ts is in O \ U, T^s e I 
and even T~^s e /, as the point is assumed not to belong to the above 
described invariant set. Moreover, by (3.14) and (3.18) 

(3.19) Ts = T(s,2) = (si,2) and T-^s ^ T~\s,2) ^ (si,l), 

thus the trajectory essentially "turns back". The time reflection sym- 
metry of (3.19) is preserved until the next fZ-hitting time. To see this 
let A; > 1 be the smallest integer (possibly 00) such that T^s e U. Up 
to time k the orbit is alternating, i.e. by induction 

T'^-'s={s',2), T'^s={s'„l), 

r-(2'-i)s = (s',1), r-2^s = (4,2) 

for some s' G S^\f/ whenever < 2/ < A;. In particular T^'^'^^s = (s", 1) 
and T-e^-^s = {s", 2) for some s" G T-'^U \ U (note that k is always 
odd). This implies by (3.18) that the preimage of T~'^''~^^s is not in /, 
alternation ceases and 

T-''s = {s'l, 2)eU cO. 
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Moreover 

T'^s = {s'l, 2) = T-''s, 

thus the point is periodic with period 2k. This way we have shown 
that any orbit which enters U twice is periodic. 

Consider any s = (s, 2) e O and the subset of generated from 
s by the reflections (3.14) and (3.15). This subset of consists of 
two rotation orbits: the first by 7 and the second by —7, which is the 
reflected image of the first one by (3.15). 

In the non-alternating rational case there are two possibilities. On 
the one hand, if both of the above mentioned two rotation orbits avoid 
U, dynamics is that of the alternating case: s is periodic of length 
that depends only on 7 (i.e. on the physical parameters). On the other 
hand, by the invariance of U with respect to (3.15) (cf. statement (2) 
of Lemma 3.4) if one of the rotation orbits hit U, the other orbit - the 
reflected image - should hit it as well. Thus if this occurs the trajectory 
of s enters U x {2} twice and becomes periodic with even length. Since 
the longest possible period occurs when we avoid U, the length of all 
possible periods is bounded above by a constant that depends only on 
7- 

In the non-alternating irrational both rotation orbits are 

dense and thus cannot avoid U, the second of the two possibilities 
sketched above occurs for all s. To complete the proof of the theorem 
we only need to show that the length of possible periods is bounded 
above. By minimality and compactness there exists a positive integer 

such that the flrst rotation-preimages of U cover S^. This implies 
that the length of any period is less than or equal to 4N. Moreover, 
N depends only on 7 and U, i.e. on the physical parameters and the 
integrals of motion. □ 



3.3. Analysis of the skew product. In this subsection we extend 
our investigation to the flbcrs, i.e. to the question how the points of 
impact for the consecutive outer bounces follow each other on the outer 
scatterer. Let us parameterize the outer circle by the arclength param- 
eter (f (the starting point is irrelevant and may be chosen arbitrarily). 
We consider the skew product 

(3.20) F : §^ X §^ ^ §^ X F{s,ip) ^ {Tos,(p + a{s)). 

Here Tq is the (projection onto of the) flrst return map onto O for T 
from Section 3.2. Note that Tq is either the flrst or the second iterate 
of T. More precisely, ii s E U, we have Tq — T and a{s) — n — 2/3(s) 
where tan/3(s) = ^ = 4^. On the other hand, if s e §^\C/, we have 

^ ^ Vn{s) W{S) ' \ ' 

To — and a{s) can be calculated using the geometry of Figure 1 
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(cf. the end of Section 2) as: 




sin/?(s) 



where 



(3.21) 



R ' 
vtjs) ^ z{s) 

Vn{s) W{s) 



tan/9(s) 



and S2 is defined by (3.15). Note that, when restricted to U or 0\U, 
both a and Tq are continuous. 

Recall from Theorem 3.7 that apart from the alternating irrational 
case To is periodic for any s with some period A'^ = ^i^)- these 
periodic situations 



F^{s,(p) = {s,(p + aN{s)); aN{s) = a{s) + a{Tos)-\ \-a{To^ ^s). 



Thus a suitable iterate of the dynamics in the fibers is a rotation, the 
angle of which depends only on the initial velocity s, and not on the 
initial position (/?. The iterates F''o+kN with < A;o < fixed and 
A; e Z, are the shifted versions of the same rotation orbit. Moreover, by 
piecewise continuity of a and Tq, and by the uniform upper bound on 
N{s), the rotation angle ci!jv(s) depends on s in a piecewise continuous 
manner as well (at discontinuity points, however, even N{s) may jump). 

Remark 3.8. Note that in the piecewise continuity statements we have 
finitely many pieces since the set U, being defined by algebraic equations 
and inequalities, has finitely many connected components. 

From this point on we restrict to the alternating irrational case. Thus 
To is an irrational rotation, and we have 



with 7 irrational and a : — > continuous. Such systems are called 
Anzai skew products after their first investigator ([A]). There is an ex- 
tensive literature on the subject, see [F],[BM], [I], [ILM] and references 
therein. Some simple issues are even treated in general monographs on 
dynamical systems like [KH]. 

In addition to the skew product map (3.22) we arc interested in 
the dynamical behavior in real time, i.e. in a suspension fiow over our 
skew product. The ceiling function for this flow, i.e. the physical (real- 
valued) flrst return time to the outer scatterer depends only on the 
initial velocity, thus it may be denoted as t(s) with r : S-*^ — > R+. We 



(3.22) 



F{s,ip) = {s + -f,cp + a{s)) 
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may calculate r(s) by Figure 2 as 

d{s) d{s2) 



r{s) 



with 



\V[S)\ \V[S2)\ 

(3.23) \v{s)\ = ^/vt{sy + v„(s)2 = y/z{sy + w{sy and 



d{l3{s)) = d{s) = cos (3{s) - \l B? - sin^ I3{s), where 

tau,a(.) = ^ = iW 

and S2 is defined by (3.15). 




Figure 2. (i(/3) = cos/3 - - sin^ (3 



We recall some basic concepts and properties related to irrational 
numbers 7 and continuous functions a : — > §^ (on details see [KH] 
and [B]). 

Definition 3.9. Given r > 1 an irrational number 7 G [0, 1] is Dio- 
phantine of type r if there exists c > such that for any p,q & 

\q1 ~ p\ > cq 

A number which satisfies the above condition for some r > 1 is called 
Diophantine. Liouville numbers 7 are those irrationals that do not 
satisfy the Diophantine condition for any r. It is not hard to see that 
the set of Diophantine numbers of type r is, as a subset of [0, 1], on the 
one hand, of the first category, and, on the other hand, of full Lebesgue 
measure, with, even more, its complement having Hausdorff dimension 
r~^. Consequently, the set of Diophantine numbers is of full measure 
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and first category, while its complement, the union of Liouville numbers 
and rational ones, has zero Hausdorff dimension. 

Definition 3.10. For any continuous map a : S-*^ — > §^ there exist 
countably many continuous functions ^4 : R — > R such that A{s + 1) — 
A{s) is an integer independent of s and a{s) can be regarded as the 
restriction of A{s) to [0, 1). The integer A{s + 1) — A{s) is independent 
of the lift and is called the degree of a{s). The map a{s) can be 
regarded as a real-valued continuous function defined on §^ iff its degree 
is zero. 

Definition 3.11. Two measurable maps a : E)^ ^ and a' : E>^ ^ 

are said to be cohomologous (for the rotation by iff there exists 
V : E)^ ^ measurable such that 

a{s) = a'{s) + V{s + ^) -V{s) 

for almost all s e . Maps cohomologous to are called coboundaries. 

The following statements date back to [A] . 

Lemma 3.12. (1) The skew product (3.22) is ergodic (with respect 
to Lebesgue measure) iff there is no p & 'L \ {0} such that pa{s) is a 
coboundary. 

(2) The skew product (3.22) has pure point spectrum iff' a{s) is not 
cohomologous to a constant. 

(3) A constant function a{s) — a is a coboundary iff a and 7 are 
rationally dependent. 

(4) By (1-3) for a function a{s) cohomologous to a constant a the 
skew product is either not ergodic or ergodic with pure point spectrum, 
depending on whether a and 7 are rationally dependent or not. 

The following result is a consequence of Theorem 1 from [BM] com- 
bined with Lemma 3.12. 

Lemma 3.13. Consider an Anzai skew product (3.22) with 7 Dio- 
phantine and a : E^ ^ E^ of zero degree and C°°. Then there exists a 
number a e [0, 1) such that the map F is smoothly (i.e. C°°) conjugate 
to 

(3.24) Fo(§^ X §1) ^ (§1 X §^), Fo{s,ip) = {s + ^,ip + a). 

We prepare for our main result related to the full system, Theorem 
3.15, with the proof of the following Lemma. 

Lemma 3.14. (1) The map a{s) is of zero degree, thus can be regarded 
asa:E^ ^R. 

(2) Both a{s) and r{s) are C°°. 

Proof. (1) Note that any continuous map S-*^ — > §^ that has nonzero 
degree must be onto. On the other hand, as we are in the alternating 
case we have cos l3{s) > R for any s e S-*^ (see Figure 1). This means 
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\c({s)\ < 2arccosi? which in particular imphes a{s) cannot be onto, 
thus it is of zero degree. 

(2) The maps a{s) and r(s) are defined in Formulas (3.21) and (3.23) 
respectively (see also (3.15)). To establish the desired smoothness prop- 
erties we need to show that four functions, namely (i) f3{s) (ii) v{s) (iii) 
/3(/9) (iv) d{l3) are all smooth. Moreover we need to show that (v) \v{s)\ 
is bounded away from zero. 

First note x, y, z all depend on s smoothly, while w a priori only 
continuously. To see this note that x{s),y{s), z{s) are the Cartesian 
coordinates of a circle (of the velocity circle, cf. Section 3.2) thus 
depend smoothly on the arclength parameter. By (3.10) 'w{s) is the 
square root of a smooth non-negative function, thus it is continuous, 
and smooth if bounded away from zero, a fact we establish below. 

By Remark 3.1 (2) and compactness 1^(5)1 is bounded away from 
zero, thus we have statement (v). This by alternating (cf. (3.17)) gives 
positive lower bound and, by the above reasoning, smoothness for w{s) 
which means (cf. (3.8)) we have shown (ii). Moreover (i) follows as we 
have demonstrated that tan f3{s) is smooth. 

We need to prove (iii) and (iv). As we are in the alternating case, 
see the first part of Remark 3.6, we have Formula (3.17) and thus by 
compactness 

sin^ 3 

- sin^ /3(s) > Ki>0; and > > 1- 

This means that the quantity under the square root in the formula for 
d{(3) (cf. (3.23)) is bounded away from zero. Similarly, in the formula 
for /9(/?) (see (3.21)), the argument of arc sin is bounded away from 1 
and —1. Thus we have (iii) and (iv). □ 



Theorem 3.15. Let us fix 7 Diophantine. 

(1) The map (3.22) is smoothly conjugate to (3.24) with a — Jgi a{s)ds. 

(2) There is a positive measure set of integrals of motion for which 
(3.22) is uniquely ergodic. 

(3) The system in "physical" time, i.e. the suspension flow over (3.22) 
with ceiling function (3.23) is never ergodic. 

Proof. (1) Lemmas 3.14 and 3.13 imply that (3.22) is smoothly con- 
jugate to (3.24). 

(2) For a and 7 rationally independent the system (3.24) is well- 
known to be uniquely ergodic, see [KH]. What we are going to demon- 
strate is that this happens for a positive measure set of parameters. 
To see this, first note that the velocity circle is determined by the 
two parameters and E (see Equation (3.9) and the discussion at 
the beginning of Section 3.2). Furthermore, it is possible to choose 
these two parameters from a positive measure set in to ensure that 
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x{s) = Vt{s) > for all s e S^. We keep the physical parameters (and 
thus 7) together with two integrals of motion, with N and E from the 
above set fixed and vary the third one, F from (3.10). What is shown 
below is that in this setting a depends on F in a continuous and strictly 
monotonic manner. This implies that a and 7 can be rationally de- 
pendent only for countably many values of F, and thus by our choice 
of N and E for a positive measure set of integrals. 

In the rest of the proof, all functions of s (e.g. z{s). f3{s),a{s)) 
refer to the original values of the integrals, while the primed ones 
(e.g. z' (s) , P' (s) , a' (s)) indicate maps for F slightly increased. Since 
Lebesgue measure on the circle is invariant with respect to the trans- 
formation s — > S2, we have 

a = [ 0{s) + P{s2))ds = 2 [ P{s)ds. 

Note that x{s) = x'{s) and y{s), z{s) similarly, however as F is 
increased we get w'{s) > w{s) for all s e S^. This implies /3'(s) < 
/3(s) (recall our convention on x{s) > 0) which yields, by a simple 
geometrical calculation, P'{s) < ${s) for all s e S-^. Thus a depends 
on F strictly monotonically. 

(3) A suspension flow above an ergodic map is ergodic iff the ceil- 
ing function is not cohomologous to a constant, see [KH]. As r(s) is 
C°° by Lemma 3.14 and 7 is chosen to be Diophantine, we may apply 
Lemma 3.13 (note also Definition 3.10) to see that r(s) is smoothly 
cohomologous to its average, a constant, thus the suspension flow is 
not ergodic. □ 

Remark 3.16. Note that an Anzai skew product (3.20) is never mixing 
as it has a factor which is a rotation by 7 (with the semi-conjugacy being 
projection to the first E)^ -coordinate). On the other hand, in case the 
spectrum is not pure point only "horizontal" sets B' = B x do not 
mix weakly, see Lemma 3.12. To be more precise, if the spectrum is not 
pure point, restricting the unitary transformation of C'^{^^ x S^) to the 
orthocomplement of the suhspace generated by characteristic functions 
of horizontal sets, the spectrum is continuous (see [A]). Even though 
this phenomenon, if occurs at all in our system, can only occur for a 
zero measure set of parameters (in the non- Diophantine case), it seems 
to he an interesting point for further investigation. 

4. Two PARTICLES, ONE ROTATING SCATTERER 

In this final section wc will show that the motion of two particles 
with one rotating scatterer, although more complicated, resembles the 
motion of one particle with two rotating scatterers. Let the two point 
particle velocity vectors be v and u. Note that as we assume that only 
the inner scatterer rotates, the outer interactions are billiard collisions 
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and thus for both particles the argument of Lemma 2.1, statement (1) 
apphes. This means that both of the particles either alternate (all odd 
bounces arc in and all even ones are out) or collide always with the 
outer billiard scatterer. We may assume that both particles alternate. 
Otherwise motion is a direct product of two integrable systems: two 
circle billiards if both particles remain outside and a product of a circle 
biUiard with the system discussed in Section 2 if one of them alternates. 

Recall Remark 1.1: in the present section splitting of v {u) into tan- 
gential and normal components {vt,Vn) {{ut,Un)) is always understood 
at the inner scatterer. Again by the geometrical picture discussed in 
the proof of Lemma 2.1 both Un and Vn are integrals of motion. For 
brevity wc drop subscripts and denote Vt {ut) by v (u). 

We need to introduce one more velocity coordinate uj for the angular 
velocity of the inner scatterer. Furthermore, as in Section 2, ij — -^^^ 
where is the moment of inertia of the scatterer and m is the mass of 
one of the point particles (we assume them to have the same mass). 

An interaction between the particle with velocity v and the inner 
scatterer will be referred to as a v-collision, and the definition of u- 
collision is analogous. In course of a v-collision the velocities transform 
as 

2r] 

v' — V (v — Ru)); 

l + rj " 

2 

(4.1) Buj' = i?a; + (v - Bw). 

We have 

(4.2) V + r]Rjj v' + riBw'\ + T]I^uy' = v'"^ + ^T^J"^ 

while u does not change. For u-collisions 

2?7 

l + ?7 
2 

(4.3) Bu)' = i?a;+ y^(m-^), 

and we have conserved quantities: 

(4.4) u + riEw = -u' + riEw'; + ryi?^^^ = -u'^ + 

in addition to the unchanged v. 

Let us introduce q = ^JfjEio. Formulas (4.2) and (4.4) imply the 
presence of two integrals: 

+ + ^ E and v + u + ^/fjq = N. 

Thus velocity motion is restricted to the intersection of a sphere and 
a plane in IR^, a circle which will be referred to as the velocity circle. 
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The V- and u-colhsions, i.e. the transformations (4.1) and (4.3) are 
reflections across the planes 

^/riv — g = and ^/r]u — q = 0, 

respectively. When restricted to the velocity circle, v-collisions reduce 
to reflections across the line 

(4.5) {N + t, -{l + v)t, V^t) 

while u-coUisions to reflections across 

( -(1 + 77)^, N + t, ^t). 

The angle of the two lines, which we denote by ^, is an important 
parameter depending on which we distinguish rational/irrational cases. 
We may calculate 

7 1 
cos — 



2 1 + 7]' 

thus, just like in Lemma 3.2, for a full measure set of rj parameters 
- actually, apart from countably many points - the base is irrational, 
and this classification does not depend on the integrals of motion. 

We parameterize the velocity circle with the normalized arclcngth 
parameter s e [0,1) as in Convention 3.3. For definiteness s = 
corresponds to one of the intersection points of the circle with the line 

(4.5) . In this coordinate the v- and u-coUisions, i.e. the transformations 
(4.1) and (4.3) are 

(4.6) s si — —s (mod 1) and 

(4.7) s S2 = -s - 7 (modi), 
respectively. 

However, at this point the exposition ceases to be parallel to that of 
Section 3.2. The two transformations above may follow each other in 
a more complicated manner and it is not enough to take two copies of 
the circle to describe the velocity evolution. 

Consider the Poincare section for the cross section of v-collisions. 
When the system is in state s, there is a unique tangential velocity 
v{s) - and thus, as the normal velocity component Vn is an integral of 
motion, a unique velocity vector v{s) - for the first particle. Thus the 
return time for the first particle can be calculated (cf. Fig 2) as 

(4.8) n{s)=T{v{s)), where t{v) = 2^ 

\v\ 

and d{v) is given in Formula (3.23) (note this time f„ does not de- 
pend on s as it is an integral of motion). The evolution of the system 
strongly depends on the number of inner bounces the second particle 
has between two consecutive bounces of the first one. As within this 
time interval the second particle follows the pattern described in Sec- 
tion 2, we have the following picture. Let us assume the system is 
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in state s e just after a v-coUision and thus the angular velocity 
of the scatterer is a;(s). If the number of u-coUisions before the next 
v-collision is even, the scatterer, just before the next v-collision, will 
rotate with the exact same angular velocity uj{s). This implies that 
just after the second v-collision the system arrives at the state Si = —s 
(cf. (4.6)). On the other hand, if the number of u-coUisions in this 
time interval is odd, the scatterer rotates with angular velocity co{s2) 
just before the second v-collision, and thus the system lands just after 
the second v-collision at the state —S2 = s + 7 (cf. (4.7)). 

To calculate the number of u-coUisions within the time interval [0, ri] 
we denote by 

(4.9) T2{s) = t{u{s)) and r^is) = t{u{s2)) 

the length of time intervals between consecutive u-coUisions in the two 
possible phases of the second particle (cf. also (4.8)). We introduce 

(4.10) A= e[0,l) and t^l \ e [0, 1) 

T2 + 7-3 t T2 Ta J 

(where {x} is the fractional part of a; G M). If the two particles started 
from the inner scatterer at the same time, the parity of u-coUisions 
until the next v-collision would depend only on the relation of these 
two quantities: it would be even/odd for i less than (greater than) A. 
To track the velocity evolution in the general case we need to introduce 
one more coordinate. 

Definition 4.1. At the moment of a v-collision (i.e. on our Poincare 
section) we define t e [0, 1) as the following quantity: the last u- 
collision preceded the v-collision we are considering at time tT2{s). 

The coordinate t defined this way describes where the second particle 
is at the moment the first particle collides inside. 

Convention 4.2. We think of the quantities t, t as elements of Ei^ 
(identified with [0, 1)) thus we may add them (which is always under- 
stood (mod 1)). 

This way the velocity evolution can be described as the following two 
dimensional system 

T : §^ X §^ ^ §^ X S\ T(s, t) = (s', t') where 

(4.11) (s', t') = ^-s, if t\ + i<\ (even case) 

(s', t') = + 7 , + if t\ + i>\ (odd case) . 

Note that T is well defined since i and A both depend on s, see Equa- 
tion (4.10). 
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Remark 4.3. The reason for the rescalinq term, in the evolution 



of the t-coordinate for the odd case is that, as there has been an odd 
number of u-collisions, the length of the period for the second particle 

is Ts in contrast to the previous T2- 

We have t' = *^^(^)+*(^2(s)+r3(s)) hoth the odd and the even cases. 

T2{S') -I 

Note u{s) = u(—s) (as the reflection (4.6) corresponds to the inter- 
action of the first particle with the scatterer, thus the velocity of the 
second particle does not change) and thus T2{s) = T2(— s) (cf (4.8) and 
(4.9) ). This implies T2{—s) — T2{s) (even case) and T2{s + 7) = 73(5) 
(odd case). 

Note that this way we need to study a two dimensional system for 
velocity evolution, which is much more complicated than the one di- 
mensional base system of Section 3.2. However, as the theorem below 
shows, for a suitable choice of integrals and, correspondingly, initial 
conditions (both of positive measure), there is an arbitrary long time 
interval for which the t coordinate can be disregarded and for which 
the velocity motion follows a rotation orbit. 

Theorem 4.4. For any e > small enough there exists a positive 
integer N^, an open set W in the space of parameters (i.e. integrals of 
motion) and a set of initial conditions E of Lebesgue measure bigger 
than s with the following property. Given a system with parameters 
from W, for any (sq^^o) £ E the iterate {sk,tk) = T^(so,to) is in 
the odd part of the phase space, and thus, Sk — Sq + k^, as long as 
< k < N^. Moreover, ^ 00 as e ^ 0. 

Proof. Let us define E C E>^ xS^ with |to~|| < £ and sq E arbitrary. 
This means that the second particle is close to the outer (billiard) 
scatterer when the first particle collides with the inner (rotating) one. 

To define the open set W let us fix and E arbitrary (in such a way 
that the velocity circle is not empty). Note that the normal velocity 
components are integrals of motion as well. We fix > small enough 
and K > big enough, their values will be chosen later depending on 
e. For the points of W we assume \vn — K'\ < e' and \un — K'\ < e' 
for some K' > K. This ensures that the normal velocity components 
of the particles are almost equal and, moreover, much greater than the 
tangential components. As a consequence, for all s, both d{v{s)) and 
d{u{s)) are very close to 1 — i?, while both \v{s)\ and \u{s)\ are very 
close to K' . In particular, with e' sufficiently small and K sufficiently 
big we may ensure that 



(4.12) 





£, 1 + e) , 



and 



t{s) 



< e 



1 - X{s) 
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for all s G S^, cf. (4.10). Note that in the last inequality the absolute 
value is understood as distance from on S^. 
Finally, we fix 

.r \ ■ ^log(2-2£) log(8e) 
(4.13) = mm ' 



log(l + 3e) ' log(l - 3e) 
where [x] denotes the integer part of a; e R. 
Lemma 4.5. Assume < k < N^. We have 

(2) X <tkX + i <1. 



[1) ^{l-3sr'<tk<l{l + 3er\ 



Note that (2) means in particular that {sk,tk) is in the odd part of the 
phase space, cf. (4.11). 



Proof. Wc will prove (1) and (2) by induction. As |to — ^| < both 
statements hold for = by (4.12). Assume that we have already 
estabhshed (1) and (2) for /c — 1. Note that this, in particular, means 
that (sfc_i,tfe_i) is in the odd part and thus by (4.11) 



tk 



^'t + ^ 

T3 1 - 



A 



This implies, by the inductive assumption and (4.12) 

ifc<^(l + 3£)^(l+£)+£ < l(l + 3£)^+\ 



tu > -(l-3£)'=(l-£) 



e > 



3e) 



fc+i 



which is (1) for tk- On the other hand, (1) implies (2) as long as k < N^. 
To see this note that by (4.12) it is enough to prove that 



-(l-3e)'=+^^(l-2£) + --£ > 



2 



2+^ 



and 



^[l + 3er^'-{l + 2e) + - + e < 1. 

These conditions are both satisfied as long as k < N^, cf. (4.13). □ 
By Formula (4.11), Theorem 4.4 is a straightforward consequence of 
Lemma 4.5. □ 



Remark 4.6. Considering motion in the fibers, i.e. in addition to the 
velocity coordinate s, the consecutive points of impact of v- collisions 
(points of impact for the first particle) parametrized by the arclength 
parameter (p along the inner scatterer, we can make the following ob- 
servation. Under the conditions of Theorem 4-4! the coordinate t can 
be disregarded, and the two dimensional system of the (s, ^p) variables 
follows (up to time N^) the Anzai skew product dynamics analyzed in 
Section 3.3. 
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